Analysis and Prediction of Urban Surface Transformation Based on Small Baseline Subset Interferometric Synthetic Aperture Radar and Sparrow Search Algorithm-Convolutional Neural Network-Long Short-Term Memory Model.

With the acceleration of urbanisation, urban areas are subject to the combined effects of the accumulation of various natural factors, such as changes in temperature leading to the thermal expansion or contraction of surface materials (rock, soil, etc.) and changes in precipitation and humidity leading to an increase in the self-weight of soil due to the infiltration of water along the cracks or pores in the ground. Therefore, the subsidence of urban areas has now become a serious geological disaster phenomenon. However, the use of traditional neural network prediction models has limitations when examining the causal relationships between time series surface deformation data and multiple influencing factors and when applying multiple influencing factors for predictive analyses. To this end, Sentinel-1A data from March 2017 to February 2023 were used as the data source in this paper, based on time series deformation data acquired using the small baseline subset interferometric synthetic aperture radar (SBAS-InSAR) technique. A sparrow search algorithm-convolutional neural network-long short-term memory (SSA-CNN-LSTM) neural network prediction model was built. The six factors of temperature, humidity, precipitation, and ground temperature at three different depths below the surface (5 cm, 10 cm, and 15 cm) were taken as the input of the model, and the surface deformation data were taken as the output of the neural network model. The correlation between the spatial and temporal evolution characteristics of the ground subsidence in urban areas and various influencing factors was analysed using grey correlation analysis, which proved that these six factors contribute to some extent to the deformation of the urban surface. The main urban area of Hohhot City, Inner Mongolia Autonomous Region, was used as the study area. In order to verify the efficacy of this neural network prediction model, the prediction effects of the multilayer perceptron (MLP), backpropagation (BP), and SSA-CNN-LSTM models were compared and analysed, with the values of the correlation coefficients of the feature points of A1, B1, and C1 being in the range of 0.92, 0.83, and 0.93, respectively. The results show that compared with the traditional MLP and BP neural network models, the SSA-CNN-LSTM model achieves a higher performance in predicting time series surface deformation data in urban areas, which provides new ideas and methods for this area of research.

A New Permanent Scatterer Selection Method Based on Gaussian Mixture Model for Micro-Deformation Monitoring Radar Images.

The micro-deformation monitoring radar is usually based on Permanent Scatterer (PS) technology to realize deformation inversion. When the region is continuously monitored for a long time, the radar image amplitude and pixel variance will change significantly with time. Therefore, it is difficult to select phase-stable scatterers by conventional amplitude deviation methods, as they can seriously affect the accuracy of deformation inversion. For different regions studied within the same scenario, using a PS selection method based on the same threshold often increases the size of the deformation error. Therefore, this paper proposes a new PS selection method based on the Gaussian Mixture Model (GMM). Firstly, PS candidates (PSCs) are selected based on the pixels' amplitude information. Then, the amplitude deviation index of each PSC is calculated, and each pixel's probability values in different Gaussian distributions are acquired through iterations. Subsequently, the cluster types of pixels with larger probability values are designated as low-amplitude deviation pixels. Finally, the coherence coefficient and phase stability of low-amplitude deviation pixels are calculated. By comparing the probability values of each of the pixels in different Gaussian distributions, the cluster type with the larger probability, such as high-coherence pixels and high-phase stability pixels, is selected and designated as the final PS. Our analysis of the measured data revealed that the proposed method not only increased the number of PSs in the group, but also improved the stability of the number of PSs between groups.

An Airborne Arc Array Synthetic Aperture Radar Vibration Error Compensation Method.

Compared to conventional radars, arc array synthetic aperture radar (SAR) enables wide-area observation under ideal conditions. However, helicopters carrying arc array SAR platforms are generally smaller in size and more sensitive to vibration, which has a greater impact on the imaging quality. In this paper, the vibration error of the arc array SAR platform is investigated, and a vibration error model of the arc array SAR platform is established. Based on the study of the vibration error model, a vibration phase estimation and compensation algorithm based on the delayed conjugate multiplication method is proposed. In the first step, distance pulse pressure processing is performed on the echo signal. In the second step, the pulse pressure signals and their delays in the same distance unit are subjected to conjugate multiplication, and the phase of the signal after conjugate multiplication is extracted. The extracted phase is then amplitude- and phase-compensated to estimate the vibration phase. In the third step, the vibration phase is compensated in the azimuthal direction of the distance pulse pressure signal, and the pairwise echo is eliminated, which completes the compensation of the airborne arc array SAR vibration platform.

Research on Time Series Monitoring of Surface Deformation in Tongliao Urban Area Based on SBAS-PS-DS-InSAR.

As urban economies flourish and populations become increasingly concentrated, urban surface deformation has emerged as a critical factor in city planning that cannot be overlooked. Surface deformation in urban areas can lead to deformations in structural supports of infrastructure such as road bases and bridges, thereby posing a serious threat to public safety and creating significant safety hazards. Consequently, research focusing on the monitoring of urban surface deformation holds paramount importance. Interferometric synthetic aperture radar (InSAR), as an important means of earth observation, has all-day, wide-range, high-precision, etc., characteristics and is widely used in the field of surface deformation monitoring. However, traditional solitary InSAR techniques are limited in their application scenarios and computational characteristics. Additionally, the manual selection of ground control points (GCPs) is fraught with errors and uncertainties. Permanent scatterers (PS) can maintain high interferometric coherence in man-made building areas, and distributed scatterers (DS) usually show moderate coherence in areas with short vegetation; the combination of DS and PS solves the problem of manually selecting GCPs during track re-flattening and regrading, which affects the monitoring results. In this paper, 45 Sentinel-1B data from 16 February 2019 to 14 December 2021 are used as the data source in the urban area of Horqin District, Tongliao City, Inner Mongolia Autonomous Region, for example. A four-threshold (coherence coefficient threshold, FaSHPS adaptive threshold, amplitude divergence index threshold, and deformation velocity interval) GCPs point screening method for PS-DS, as well as a Small Baseline Subset-Permanent Scatterers-Distributed Scatterers-Interferometric Synthetic Aperture Radar (SBAS-PS-DS-InSAR) method for selecting PS and DS points as ground control points for orbit refinement and re-flattening, are proposed. The surface deformation results obtained using the Small Baseline Subset Interferometric Synthetic Aperture Radar (SBAS-InSAR) and the SBAS-PS-DS-InSAR proposed in this paper were comparatively analysed and verified. The maximum cumulative line-of-sight settlements were -90.78 mm and -83.68 mm, and the maximum cumulative uplifts are 74.94 mm and 97.56 mm, respectively; the maximum annual average line-of-sight settlements are -35.38 mm/y and -30.38 mm/y, and the maximum annual average uplifts are 25.27 mm/y and 27.92 mm/y. The results were evaluated and analysed in terms of correlation, mean absolute error (MAE), and root mean square error (RMSE). The deformation results of the two InSAR methods were evaluated and analysed in terms of correlation, MAE, and RMSE. The errors show that the Pearson correlation coefficients between the vertical settlement results obtained using the SBAS-PS-DS-InSAR method and the GPS monitoring results were closer to 1. The maximum MAE and RMSE were 13.7625 mm and 14.8004 mm, respectively, which are within the acceptable range; this confirms that the monitoring results of the SBAS-PS-DS-InSAR method were better than those of the original SBAS-InSAR method. SBAS-InSAR method, which is valid and reliable. The results show that the surface deformation results obtained using the SBAS-InSAR, SBAS-PS-DS-InSAR, and GPS methods have basically the same settlement locations, extents, distributions, and temporal and spatial settlement patterns. The deformation results obtained using these two InSAR methods correlate well with the GPS monitoring results, and the MAE and RMSE are within acceptable limits. By comparing the deformation information obtained using multiple methods, the surface deformation in urban areas can be better monitored and analysed, and it can also provide scientific references for urban municipal planning and disaster warning.

System Design and Echo Preprocessing of Spaceborne Squinted Two-Dimensional Beam Scanning Synthetic Aperture Radar.

Conventional squinted sliding spotlight synthetic aperture radar (SAR) imaging suffers from substantial swath width reduction and complex processing requirements due to the continuous variation in the squint angle and the large range cell migration (RCM) throughout the data acquisition interval. A novel two-dimensional (2D) beam scanning mode for high-resolution wide swath (HRWS) imaging is proposed. The key to the novel imaging mode lies in the synchronous scanning of azimuth and range beams, allowing for a broader and more flexible imaging swath with a high geometric resolution. Azimuth beam scanning from fore to aft was used to improve the azimuth resolution, while range beam scanning was adopted to illuminate the oblique wide swath to avoid the large RCM and the serious swath width reduction. Compared with the conventional sliding spotlight mode, both the swath width and swath length could be extended. According to the echo model of this imaging mode, an echo signal preprocessing approach is proposed. The key points of this approach are range data extension and azimuth data upsampling. A designed system example with a resolution of 0.5 m, swath width of 60 km, and azimuth coverage length of 134 km is presented. Furthermore, a simulation experiment on point targets was carried out. Both the presented system example and imaging results of point targets validated the proposed imaging mode.

A Modified Keystone-Based Forward-Looking Arc Array Synthetic Aperture Radar 3D Imaging Method.

An arc array synthetic aperture radar (AA-SAR) is a new type of omnidirectional observation and imaging system. Based on linear array 3D imaging, this paper introduces a keystone algorithm combined with the arc array SAR 2D imaging method and proposes a modified 3D imaging algorithm based on keystone transformation. The first step is to discuss the target azimuth angle, retain the far-field approximation method of the first-order term, analyze the influence of the forward motion of the platform on the along-track position, and realize the two-dimensional focusing of the target slant range-azimuth direction. The second step is to redefine a new azimuth angle variable in the slant-range along-track imaging and use the keystone-based processing algorithm in the range frequency domain to eliminate the coupling term generated by the array angle and the slant-range time. The corrected data are used to perform along-track pulse compression to obtain the focused image of the target and realize the three-dimensional imaging of the target. Finally, in this article, the spatial resolution of the AA-SAR system in the forward-looking state is analyzed in detail, and the change in the spatial resolution of the system and the effectiveness of the algorithm are verified through simulation.

Azimuth Full-Aperture Processing of Spaceborne Squint SAR Data with Block Varying PRF.

The block varying pulse repetition frequency (BV-PRF) scheme applied to spaceborne squint sliding-spotlight synthetic aperture radar (SAR) can resolve large-range cell migration (RCM) and reduce azimuth signal non-uniformity. However, in the BV-PRF scheme, different raw data blocks have different PRFs, and the raw data in each block are insufficiently sampled. To resolve the two problems, a novel azimuth full-aperture pre-processing method is proposed to handle the SAR raw data formed by the BV-PRF scheme. The key point of the approach is the resampling of block data with different PRFs and the continuous splicing of azimuth data. The method mainly consists of four parts: de-skewing, resampling, azimuth continuous combination, and Doppler history recovery. After de-skewing, the raw data with different PRFs can be resampled individually to obtain a uniform azimuth sampling interval, and an appropriate azimuth time shift is introduced to ensure the continuous combination of the azimuth signal. Consequently, the resulting raw data are sufficiently and uniformly sampled in azimuth, which could be well handled by classical SAR-focusing algorithms. Simulation results on point targets validate the proposed azimuth pre-processing approach. Furthermore, compared with methods to process SAR data with continuous PRF, the proposed method is more effective.

An Adjusted Frequency-Domain Algorithm for Arc Array Bistatic SAR Data with One-Moving Transmitter.

Arc array synthetic aperture radar (AA-SAR), which can observe the scene in all directions, breaks through the single view of traditional SAR. However, the concealment of AA-SAR is poor. To mitigate this, arc array bistatic SAR (AA-BiSAR) with the moving transmitter is proposed, it has the advantages of good concealment and can expand the imaging scene, and improve the flexibility of the system. The imaging geometry including the signal model is established, and a range frequency-domain algorithm based on keystone transform (KT) is proposed in this paper. In the first step, the slant range equation is approximated by Taylor series expansion to compensate for the residual phase caused by the transmitter motion. In the second step, the range cell migration between the range and azimuth is eliminated through the KT method in the range frequency-domain. In the third step, using the data after range cell migration correction in step 2, an azimuth pulse compression is performed to obtain the focused image. In addition, the spatial resolution of the AA-BiSAR system is analyzed in detail. Finally, three simulation results verify the effectiveness of the proposed algorithm and the change in the spatial resolution.

Staring Spotlight SAR with Nonlinear Frequency Modulation Signal and Azimuth Non-Uniform Sampling for Low Sidelobe Imaging.

In synthetic aperture radar (SAR) imaging, geometric resolution, sidelobe level (SLL) and signal-to-noise ratio (SNR) are the most important parameters for measuring the SAR image quality. The staring spotlight mode continuously transmits signals to a fixed area by steering the azimuth beam to acquire azimuth high geometric resolution, and its two-dimensional (2D) impulse response with the low SLL is usually obtained from the 2D weighted power spectral density (PSD) by the selected weighting window function. However, this results in the SNR reduction due to 2D amplitude window weighting. In this paper, the staring spotlight SAR with nonlinear frequency modulation (NLFM) signal and azimuth non-uniform sampling (ANUS) is proposed to obtain high geometric resolution SAR images with the low SLL and almost without any SNR reduction. The NLFM signal obtains non-equal interval frequency sampling points under uniform time sampling by adjusting the instantaneous chirp rate. Its corresponding PSD is similar to the weighting window function, and its pulse compression result without amplitude window weighting has low sidelobes. To obtain a similar Doppler frequency distribution for low sidelobe imaging in azimuth, the received SAR echoes are designed to be non-uniformly sampled in azimuth, in which the sampling sequence is dense in middle and sparse in both ends, and azimuth compression result with window weighting would also have low sidelobes. According to the echo model of the proposed imaging mode, both the back projection algorithm (BPA) and range migration algorithm (RMA) are modified and presented to handle the raw data of the proposed imaging mode. Both imaging results on simulated targets and experimental real SAR data processing results of a ground-based radar validate the proposed low sidelobe imaging mode.

Time- and Space-Varying Atmospheric Phase Correction in Discontinuous Ground-Based Synthetic Aperture Radar Deformation Monitoring.

Ground-based synthetic aperture radar (GB-SAR) uses active microwave remote-sensing observation mode to achieve two-dimensional deformation measurement and deformation trend extraction, which shows great prospects in the field of deformation monitoring. However, in the process of GB-SAR deformation monitoring, the disturbances caused by atmospheric effect cannot be neglected, and the atmospheric phases will seriously affect the precision of deformation monitoring. In discontinuous GB-SAR deformation monitoring mode, the atmospheric phases are particularly affected by changes of time and space, so the traditional models of atmospheric phase correction are no longer applicable. In this paper, the interferometric phase signal model considering atmospheric phase is first established. Then, the time- and space-varying characteristics of the atmospheric phase are analyzed, and a novel time- and space-varying atmospheric phase correction algorithm, based on coherent scatterers analysis, is proposed. Finally, slope deformation monitoring experiments are carried out to verify the validity and robustness of the proposed algorithm.

δ-Generalized Labeled Multi-Bernoulli Filter Using Amplitude Information of Neighboring Cells.

The amplitude information (AI) of echoed signals plays an important role in radar target detection and tracking. A lot of research shows that the introduction of AI enables the tracking algorithm to distinguish targets from clutter better and then improves the performance of data association. The current AI-aided tracking algorithms only consider the signal amplitude in the range-azimuth cell where measurement exists. However, since radar echoes always contain backscattered signals from multiple cells, the useful information of neighboring cells would be lost if directly applying those existing methods. In order to solve this issue, a new δ-generalized labeled multi-Bernoulli (δ-GLMB) filter is proposed. It exploits the AI of radar echoes from neighboring cells to construct a united amplitude likelihood ratio, and then plugs it into the update process and the measurement-track assignment cost matrix of the δ-GLMB filter. Simulation results show that the proposed approach has better performance in target's state and number estimation than that of the δ-GLMB only using single-cell AI in low signal-to-clutter-ratio (SCR) environment.

Multichannel High Resolution Wide Swath SAR Imaging for Hypersonic Air Vehicle with Curved Trajectory.

Synthetic aperture radar (SAR) equipped on the hypersonic air vehicle in near space has many advantages over the conventional airborne SAR. However, its high-speed maneuvering characteristics with curved trajectory result in serious range migration, and exacerbate the contradiction between the high resolution and wide swath. To solve this problem, this paper establishes the imaging geometrical model matched with the flight trajectory of the hypersonic platform and the multichannel azimuth sampling model based on the displaced phase center antenna (DPCA) technology. Furthermore, based on the multichannel signal reconstruction theory, a more efficient spectrum reconstruction model using discrete Fourier transform is proposed to obtain the azimuth uniform sampling data. Due to the high complexity of the slant range model, it is difficult to deduce the processing algorithm for SAR imaging. Thus, an approximate range model is derived based on the minimax criterion, and the optimal second-order approximate coefficients of cosine function are obtained using the two-population coevolutionary algorithm. On this basis, aiming at the problem that the traditional Omega-K algorithm cannot compensate the residual phase with the difficulty of Stolt mapping along the range frequency axis, this paper proposes an Exact Transfer Function (ETF) algorithm for SAR imaging, and presents a method of range division to achieve wide swath imaging. Simulation results verify the effectiveness of the ETF imaging algorithm.

Investigation of Azimuth Multichannel Reconstruction for Moving Targets in High Resolution Wide Swath SAR.

The azimuth multichannel imaging scheme with the large receive antenna divided into multiple sub-apertures usually leads to azimuth non-uniform sampling, and echoes from all azimuth channels should be reconstructed based on the signal model before conventional SAR imaging. Unfortunately, the multichannel signal model of a moving target is different from that of a fixed target. This paper analyzes the multichannel signal model of the moving target and the effect of the target velocity on azimuth multichannel reconstruction. Based on the multichannel signal mode of the moving target, a new multichannel signal reconstruction algorithm is proposed. Furthermore, the slant range velocity is estimated by computing signal energy distribution. Simulation results on point targets validate the proposed approach.